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Abstract 

This paper presents a simple model that mimics quantum mechanics (QM) 
results in terms of probability fields of free particles subject to self-interference, 
without using Schrodinger equation or complex wavefunctions. Unlike the 
standard QM picture, the proposed model only uses integer-valued quanti¬ 
ties and arithmetic operations. In particular, it assumes a discrete spacetime 
under the form of an euclidean lattice. The proposed approach describes in¬ 
dividual particle trajectories as random walks. Transition probabilities are 
simple functions of a few quantities that are either randomly associated to 
the particles during their preparation, or stored in the lattice sites they 
visit during the walk. Non-relativistic QM predictions, particularly self¬ 
interference, are retrieved as probability distributions of similarly-prepared 
ensembles of particles. Extension to interacting particles is discussed but 
not detailed in this paper. 

1 Introduction 

Despite being an extremely succesful theory to predict the outcome of ex¬ 
periments with particles and other microscopic objects, quantum mechanics 
(QM) is believed by many to be incomplete or, at least, not fully understood. 
In particular, the “strange” or non-classical phenomena of QM, like self¬ 
interference and Born’s rule, are described in terms of abstract mathemati¬ 
cal objects. Although some have been ready to interpret the complex-valued 
wavefunction as a real object, wavefunctions are generally seen as mathe¬ 
matical tools serving to calculate probabilities from their square moduli. 
Contrasting to real-valued mathematics and one-to-one mapping between 
real variables and observables of classical theories, the standard description 
is thus sometimes considered as purely operational. 

In this paper a model is presented that, although very simple in terms 
of mathematical development and formalism, seems to be able to predict 
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probability fields of at least free particles, and particularly self-interference, 
without using Schrodinger equation, wavefunctions, or path integrals. The 
proposed model assumes a discrete spacetime under the form of a euclidean 
lattice. Individual particle trajectories are described as random walks. Tran¬ 
sition probabilities are simple functions of a few quantities that are either 
randomly associated to the particles during their preparation, or stored in 
the lattice sites they visit during the walk. Non-relativistic QM predictions 
are retrieved as probability distributions of similarly-prepared ensembles of 
particles. Unlike the standard QM picture |3], the proposed model only 
uses integer-valued quantities and arithmetic operations, avoiding complex¬ 
valued operations on the random walks. 

The proposed model is different from ensemble interpretations mm in 
the sense that it also describes the behavior of individual particles. With 
respect to De Broglie-Bohm mechanics p] and deterministic trajectory rep¬ 
resentation [2], the proposed model introduces a non-deterministic behavior 
but does not appeal to non-localities. 

Intrinsic indeterminism is already contained in stochastic interpretations 
of QM that are, however, mostly aimed at retrieving the Schrodinger equa¬ 
tion from a classical equation of motion plus a stochastic force. The Born 
probability rule remains unexplained in this context mmm, or is founded 
on the definition of probability density of particles as the squared intensity 
of an associated wave [9|. This latter assumption is not used in the proposed 
approach, which in contrast predicts nonclassical consequences of Born rule 
(as double-slit interference) only from the random walk features. However, 
unlike other corpuscolar models m the proposed model is not restricted to 
double-slit intereference as it is not explicitly based on ad-hoc trigonometric 
functions inspired by the complex waveforms structure. 

The idea of lattice or discrete-time algorithms that reproduce particle 
propagation in the continuum limit is also not new. However, the proposed 
model uses the lattice only as the support for particle motion, not for wave- 
functions or other mathematical operators as, e.g., in PHHatH]. While 
other random walks or spacetime quantizations mm [16] are able to re¬ 
produce the emergence of Schrodinger equation from pure combinatorics, 
again the Born probability rule and thus interference are not explained in 
such models, while naturally emerges in the proposed one. In order to re¬ 
produce interference, antiparticles are not appealed to, as in some abstract 
lattice gas models for waves m, nor negative probabilities. 

The paper is organized as follows. Section 2 presents the assumptions 
concerning spacetime, which naturally lead to Heisemberg’s uncertainty 
principle. Section 3 describes the model for free motion without quantum 
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forces and retrieves de Broglie relation and Schrodinger equation. Section 4 
describes the model with quantum forces and shows numerical results for 
several self-interference scenarios. Section 5 is just a short introduction to 
possible extensions to treat interacting particles and the transition to rela¬ 
tivistic QM. 

2 The lattice: uncertainty principle 

The proposed model assumes that the spacetime is inherently discrete. Lim¬ 
iting for simplicity the analysis to one dimension x, that means that only 
values x = £X, £ € Z and t = rT, t£N are meaningful. Noninteger values 
of space and time are simply impossible in this picture. The two fundamen¬ 
tal quantities X and T are the size of the lattice that constitutes the space 
and the fundamental temporal resolution, respectively. 

Under this assumption, a particle’s history consists of a succession of 
points {£ n ,T n } in the spacetime, where n € Z is the discrete index that 
describes advance in history, here denoted as “iteration”. Advance in time 
is unidirectional and unitary, that is, r + 1 follows necessarily r. Advance 
in space is still unitary but bidirectional. If in a certain iteration a particle 
resides at the location £ of the spatial lattice, in the next iteration the 
particle can only reside at locations £ + 1, £, or £ — 1. The local velocity 
of this mechanism, V, is a random variable that can take only the three 
discrete values v = {+1, 0, —1}, as described in Sect. [3] and Sect. 01 

The maximum velocity that a particle can reach is the speed of light c. 
Light trajectory in the positive direction corresponds to v n = 1, Vn. Conse¬ 
quently, one constraint to the fundamental lattice quantities is necessarily 

X 

T 

The average velocity could only be estimated by an observer as a sample 
mean of successive v’s. Every observation lasting N iterations will give an 
approximation v of the “true” average velocity. Consider, e.g., N = 1. The 
observed velocity can be +c, 0 or —c. Thus the uncertainty on v is c in 
absolute value. For N = 2, the possible outcomes for the sample mean are 
c, c/2, 0, —c/2, and —c. Thus the uncertainty on v is c/2 in absolute value. 
Extending these considerations, the uncertainty on v after an observation 
lasting N iterations is c/N in lattice units. 

Moreover, an observation lasting N iterations necessarily implies a change 
in the position of the particle. The span of the particle during the observa- 
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tion ranges from NX to —NX. Thus the uncertainty on the position of the 
particle at the end of the observation is 2N in lattice units. 

Using the two results above, and denoting A v(N) and A x(N) the un¬ 
certainties of velocity and position as a function of observation horizon N, 
the relationship 


A v(N) ■ A x(N) 


N 


2 NX = 


2X 2 

~T~ 


( 2 ) 


holds. 

The latter equation resembles the Heisenberg uncertainty principle since 
it fixes an inverse proportionality between the uncertainty with which the 
velocity of a particle can be known and the uncertainty with which its posi¬ 
tion can be known. Multiplying by the particle mass m. and comparing ((2|) 
to Heisenberg uncertainty principle, one obtains that the two fundamental 
lattice quantities are related to the Planck constant, 


X 2 _h _ h 

m ~T ~ 2 ~ 4tt' 


( 3 ) 


The term An r (the solid angle of a sphere) holds for three-dimensional spaces. 
In our example case of a one-dimensional space, this term reduces to 2, the 
measure of the unit 1-sphere. Thus, combining (JX|) with the accordingly 
modified ([3]), the values for the fundamental lattice quantities are obtained 
as 


and 


X = 


T = 


h 


2 me 


h 


( 4 ) 


( 5 ) 


2 me 2 

Note that the Compton wavelength is retrieved as twice the fundamental 
lattice size X. 

The role of mass is not completely clear at this point. Likely, general 
relativity will serve to integrate it into the picture. 


3 Free motion without interference 

This section will first describe the equations of motion of a free particle 
in the proposed model. Then, the stochastic variables associated with the 
motion will be analyzed. Finally, the equivalence with the wavefunction 
picture and Schrodinger equation will be retrieved. 
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3.1 Particle dynamics 

This section describes the propagation rules of a particle on the lattice, or 
its dynamics. 

3.1.1 Equations of motion 

As stated in the previous section, time can only increase by one unit at each 
iteration. Time is re-initialized to zero whenever the particle interacts with 
the environment (external forces). This event is called preparation in the 
following. A stochastic variable representing time at an iteration n could be 
introduced as 

n 

Tn ■= E L ( 6 ) 

n'=no 

where no is the most recent iteration when the particle has undergone prepa¬ 
ration, that is, has interacted with the environment and has been actual¬ 
ized. However, T n would have a deterministic distribution, Pr(T n = r) = 
5(t — (n — no)) and thus it will be often replaced by its support r = n — no 
in the following. 

Consider now spatial dynamics. At each iteration, the particle might 
jump to one of the nearest neighboring sites of the lattice, or stay at rest. 
The actual local trajectory is not deterministic, i.e., it is not a prescribed 
function of previous parts of trajectory. Rather, the local trajectory has the 
characteristics of a random walk. This point is very important and it implies 
that an intrinsic randomness affects the particle motion. Generally, there is 
a different transition probability for each of the three possible transitions. 
In free motion without interference or external forces, these probabilities 
do not change with time. Let me denote the transition probabilities a := 
Pr(V n = 1), b := Pr(V„ = 0), and c := Pr(V n = —1), respectively. Of course, 

a + b + c = 1. (7) 

Moreover, the proposed model assumes that the expected value of V„ is 
imprinted to the particle. This imprint is to be attributed to the preparation 
and is actualized every time the particle interacts with the environment (in a 
way to be considered later). Let me denote this expected value as momentum 
propensity p, 

p := E[V n ] = a - c. (8) 

Another characteristic of the random motion is the expected value of the 
squared velocity, that is, 

e := E[Vn] = a + c (9) 
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that can be reinterpreted as an energy propensity. Combining (J7 |) -l]9 |) . obtain 


a = 


e+p 

2 


b = 1 — e, 


c = 


e — p 
2 


( 10 ) 


The energy e must be a function of p. A well-known result of special 
relativity states that energy of a particle is the sum of the rest energy and 
the kinetic energy. Following this suggestion, the proposed model assumes 
that 

e(p) := l -±f. (11) 

Equation m might be also interpreted in the following way: energy is 
the average of the “time energy” and the “space energy”, where the for¬ 
mer contribution is always one, since time can only advance by one unit. 
Consequently, m can be rewritten as 


a = 




c = 


¥0 


2 


( 12 ) 


Appendix IA.1I shows that it is possible to retrieve the de Broglie relation 
E = hu: for “matter waves” with the proposed model, involving the energy 
propensity e. 


3.1.2 Probability mass functions 

The equations in the previous section, in fixing the probability of each jump 
at each iteration n, define the trajectory of the particle as a random walk. 
Let me introduce the stochastic variables that describe the position of the 
particle, X ni 

n 

Xn := 6 + Vn'- (13) 

n'=riQ 

where £o denotes the site occupied at n = no, i.e., r = 0. 

Let me calculate the probablity mass function of this stochastic variable, 

Pn(0 ■= Pr(Xn = 0, (14) 

i.e., the probability that the particle has crossed £ sites after n iterations. 
Consider the scenario where particles are emitted from a source located at 
the site £o = 0 of the lattice with an intrinsic value of p, and thus of e, 
determined by the preparation. Time interval between two emissions is very 
large, so to exclude any interactions between successive particles. Moreover, 
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the single source excludes quantum interference. After one iteration or, 
equivalently, time step, the particle has a probability a to be at the site 
£ = 1, a probability 6 to be at the site £ = 0, and a probability c to be 
at the site £ = —1. After two iterations, the probabilities are: p2(2) = a 2 , 
P 2 (1) = 2 ab, P 2 (0) = 2 ac + b 2 , 1) = 26c, /C> 2 (—2) = c 2 . Note that, 

since the functions a(p) and c(p) are symmetric, the probability function is 
symmetric with respect to £ = 0. 

In general, the position probability function is described by the recursive 
equation 

Pn(£) = apn-l(£ - 1) + 6p n _l(£) + cp n - 1 (£ + 1). (15) 

Deriving a closed formula for the probability mass function p n (£) is te¬ 
dious but straightforward at this point. With no = 0 and the initial condi¬ 
tion p o(0) = 1, the result is 

( 2T ) 

Pn(i) = -- (1 +J>) T+£ (1 (16) 

with support £ = {—r,... , r}, and can be verified by inspection. From this 
formula, the probability that a particle is at the event horizon, i.e., p n ( t), 
is easily retrieved as a T . Similarly, p n (— r) = c T . The estension to the 
more general case of £o / 0 is easily done by replacing £ with £ — £o in the 
right-hand side of (fTBll . Finally, it turns out that E[X n \ = pr + £q. 

The function (11611 has a limit for large r’s that can be derived in two 
equivalent ways. On the one hand, the equation of motion C3D can be 
reviewed in the continuum limit as a stochastic differential equation reading 

dX n = E[V n ]dn + VVar [V n ]dB, (17) 

where dB is a Brownian motion with zero mean and unit variance. Now, 
from (j8] l - ([T2l) . the identity Var[V n ] = e — p 2 = 6 follows. Consequently, the 
continuum limit of p is a Gaussian function with mean pr and variance 6r, 
that is, 

:= to ft.K) » ■ < 18 ) 

Other derivations of (11811 are illustrated in Appendix IA.2I Appendix IA.3I 
shows that (11811 is invariant to Lorentz transformations. 

Let me define an additional stochastic variable for later use, representing 
the distance covered by the particle 

n 

Cn := J2 Vn ' = Xn-Zo. (19) 

n'=riQ 
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Clearly, Pr (C n = A) = Pr(,T n = £o + A) and E[C n \ = pr hold. 


3.2 Lattice dynamics 


To describe particle dynamics as “seen” by the lattice sites, let me define 
additional stochastic variables. Denote these variables with a time subscript 
and a position superscript, instead of the only time subscript as for the 
particle variables. Clearly, the variable can only take the value £ in the 
single-source scenario under consideration. Similarly, a variable Tr would 
be deterministic, with its support including only the value r. 

Let me introduce a true stochastic variable, the site occupancy, defined 


as 


ot 


1, if the particle occupies the site {£,r} 
0, otherwise 


Define pi := Pr (Of = 1) = E[Oi\. Clearly, 


( 20 ) 


pi =pao 


( 21 ) 


and its pmf is given by (USD- 


3.3 Ensembles of Particles 

In this section the predictions of QM are retrieved for ensembles of similarly 
prepared particles, that is, particles emitted from the same source at £o = 0. 


3.3.1 Probability density 

To retrieve the predictions of Schrodinger’s equation, a key element of the 
model is introduced. The proposed model assumes that whenever the parti¬ 
cle interacts with the environment, its momentum propensity p is properly 
reset. Now, in the free motion scenario, consider for the moment that p 
is a continous variable determined randomly during the preparation at the 
particle source. Consequently, the probability of releasing a particle with 
a momentum propensity p is uniform over the interval between -1 and +1, 
spanning 2, and thus the probability density of the momentum propensity 
is f(p) = 1/2. 

When the source releases a large number of particles in succession, each 
one with a randomly determined value of p, the probability of finding a par¬ 
ticle at the location {£, r} is given by the ensemble average and is calculated 
as ^ 

Pi ■= J i f(p)pldp =\ j i PUP- ( 22 ) 
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Introducing (fl6l) into (l22l) , and after some manipulations (see Appendix lA.51) , 
obtain 

p r = ^ (23) 

Moreover, it is easily verified that 

£ P$ = 1 ( 24 ) 

as obviously required. 

Now, compare this result with the predictions of QM, i.e., the particular 
solution of the Schrodinger equation. For a single perfectly localized source 
at x = 0, the probability density is calculated fsee IA.6I) as 


l*0M)l 2 


mX 2 

2irht 


mX 2 

ht 


(25) 


thus it is inversely proportional to time and it does not depend on x. Nor¬ 
malizing to lattice units and using Q allows reducing (1251) to 


l<t(f ' T)|2= 2? (26) 

This result compares with (1231) . with 2r replacing 2r + 1. The two functions 
of r are very similar and, indeed, practically coincident for r sufficiently 
large. In other terms, the square modulus of the wavefunction predicted by 
the Schrodinger equation is the continuum limit of the probability in the 
proposed model. 

Note that only the particular formulation of energy propensity m yields 
this result. Other values for e and consequently for b (for instance, 6 = 0), 
would yield position-dependent mass probability functions. 


3.3.2 Phase 

To retrieve the phase of the wavefunction predicted by the Schrodinger equa¬ 
tion, define first a new lattice stochastic variable St representing the accu¬ 
mulated energy of the particle as seen by a site. Its pmf is calculated (see 
Appendix IA.4P as 


4 ■= Pr (Sf = a) 




(27) 
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with support a = {|£|, |£| + 2,..., |£| + 2 }. It can be also proved that 



(28) 


and 



(29) 


Equation (051) can be easily verified by inspection of a few sites, provided 
that {£,r} > 0. For example, for £o = 0, one easily obtains £{5}] = 1. 
The cumulated energy of every particle reaching the site {1,1} is obviously 
1. Generalizing this result, clearly E[S^\ = r. For a site like {0,2} the 
prediction is less trivial. The possible values of the cumulated energy can 
be 2, if the particle follows a back-and-forth path, or 0, if it stays at rest 
for two time steps. Using (|28[) . one obtains for this case FfiSf, 1 ] = 2/3, which 
is a weighted mean between the two possible values of cumulated energy. 
Equation ()29D also can be verified by inspection. For instance, when £ = ±t 
or £ = ±(t — 1), it is apparent that <5>| can only take the value |£|. In fact, 
(1291) yields Var[«S|] = 0. 

The continuum limit counterpart of (l27l) is calculated as 



(30) 


Define now the action at any site as 



(31) 


Note that the expected value given by equation (f28l) does not depend on p. 
Therefore, by inserting it into (1311) . obtain 



(32) 


Compare this result with the phase of the wavefunction (1901) . The latter, 
usually interpreted as the action of the particle, is 



(33) 
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which, in lattice units, becomes 



(34) 


The latter equation corresponds to the second term in the right-hand side 
of (1321) . that is, X$ — X^ 0 , multiplied by n to obtain a phase angle. The 
correspondence is almost perfect, except for the term 2r — 1 that in the 
proposed model replaces the term 2r predicted by QM. For large values of 
r, however, the two results are practically coincident. In other terms, the 
phase of the wavefunction predicted by the complex Schrodinger equation is 
the continuum limit of the action in the proposed model. 

3.3.3 Schrodinger equation: de Broglie-Bohm formulation 

The results in the previous sections have been derived for a probability 
density of the momentum propensity f(p) = 1/2. Consider now a generic 
function f{jp). Apply (l22l) to the continuum limit (fT8l) of the probability 
function p%. It can be shown that (1181) approximates a Dirac delta function, 
whence 



(35) 


The probability density function P obeys the following partial differential 
equation 



Introducing now 


the continuum-limit approximation of X$, 



(37) 


one recognizes in (l36l) the continuity equation 



(38) 


On the other hand, the relationship 



(39) 


also holds. 
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Sect. 13.3.11 has shown the equivalence of P to |\H| 2 , while Sect. 13.3.21 
that of X with S/ir = ZT/vr for large r’s. With these two substitutions, 
and reintroducing physical units instead of lattice units, equations (l38li - (1391) 
become the continuity equation 

__ ( m 2 dSH/rn\ 

dt dx \ dx ) 

and the Hamilton-Jacobi equation 

dS_ _ _ 1_ (dSh\ 2 

dt 2m \ dx ) ’ 

of the de Broglie-Bohm formulation of QM, which in turn are equivalent to 
Schrodinger’s equation for a free particle. 

3.4 Numerical results 

In the last sections, the predictions of the proposed model were shown in 
closed form, using mathematical equations in terms of a priori probabilities 
and probability fluxes. The probability of a number of observable were 
calculated and found to be in accord to the predictions of QM. Now, I will 
present numerical simulations of the random walk of single particles and I 
will calculate the a posteriori probabilities as frequencies over a large number 
of emissions. Thus, this section is aimed at reproducing numerically a true 
experiment. 

Table |T] shows the algorithm used for such simulations. The two for- 
cycles are for the successively released Np particles, and for time up to Np. 
Each particle experiences the choice of two randomly-selected values: (i) the 
momentum propensity and (ii) at each time step, its local velocity v as a 
function of p. The final code line represents the counting of the particle that 
arrive at a certain location at time Np. From this number of arrivals, an a 
posteriori frequency z/(£) is calculated as the ratio to the total number of 
particles emitted. 

Figure Q] shows the frequency i/(£) after a time Np = 300 for different 
values of Np. As the the number of particles emitted in the ensemble in¬ 
creases, a frequency distribution builds up. For large Np, the frequency 
clearly tends to the a priori probability P(£,r = Np), that is, a constant 
value given by (1231) . 
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Algorithm 1 Pseudocode used for the simulations of Fig. [I] 

1: for particle 1 to Np do 
2 : £<- 0 

3: T •$— 0 

4: p <r- random value beween -1 and +1 

5: for iteration 1 to Np do 

6: r 4— t + 1 

7: v <— random value +1, 0 or -1 with prob. given by (1101) 

8 : £ <T- £ + V 

9: end for 

10 : 1/(0 «- + 1 

11: end for 
12: z/(£) = v{£)/Np 



-200 0 200 



£ 


x 10" 3 



x 10^ 3 
6 -- 



2 


0 1 -*-*-— 

-200 0 200 

£ 


Figure 1: Frequency of arrival of particles emitted at £ = 0 as a function of 
£ after Nt = 300. From top-left to bottom-right, Np = 500, 5000, 10000, 
and 50000, respectively. 
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4 Interference 


After having reproduced the predictions of the Schrodinger equation for 
a free particle, let me proceed now to a second puzzling aspect of QM: 
particle self-interference. Double-slit experiment usually serves to visualize 
this phenomenon. However, the core of self-interference is isolated and better 
illustrated by a double-source preparation, where particles can be emitted by 
two alternative sources and the two possible paths interefere with each other. 
Additional free-particle scenarios leading to self-intereference are multiple- 
source preparations and the “particle in a ring” situation. 

The representation of these scenarios using only the process illustrated 
in Sect. [3] would give just the superimposition of probability densities of 
the type (1231) and no interference would arise. Instead, interference is orig¬ 
inated in the proposed model from the interaction of successive particles 
emitted with the lattice, a process denoted here as “quantum force,” which 
is illustrated in the next sections. 

4.1 Quantum force mechanism 

In short, the emergence of quantum forces is the result of an exchange of 
information between the lattice and the particles. Let me call this piece of 
information boson , in analogy with force-mediating particles. 

In the proposed model, each particle carries two counters. The spatial 
counter records the number of space steps crossed and thus it is represented 
by the stochastic variable C n . The temporal counter records the number of 
time steps crossed and thus is always certainly equal to r. 

On the other hand, each lattice site maintains a sort of “register” storing 
the value of the spatial counter of the last particle that has visited the site. 
Let me describe this register by the lattice stochastic variable C\. 

When a particle visits a site, the site senses the value of the particle’s 
spatial counter, say, A. This value is compared with the value /r stored in the 
register, which is the value sensed the last time a particle visited the site. 
If the difference 5 := \n — A| is strictly positive, an exchange of information 
between the particle and the lattice site takes place. 

First, a couple of “bosons” is created. One boson is taken by the particle 
(later it will be explained how this boson acts on the particle), while the 
second stays at the site. The pair {A, ^} for which these bosons have been 
created distinguishes them from other possible bosons. For this reason let 
me call them A/U-bosons. 

The second event taking place when 6 > 0 is that the value fi in the site 
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register is replaced with the value A, while the particle counter is replaced 
with the value /i. In other words, the site register and the particle counter 
are exchanged, 

£$ 4— A, C n <— fi. (42) 

The next subsections will describe separately the behavior of site bosons 
and particle bosons. 


4.1.1 Quantum force at the lattice 

When a new A/U-boson is created at a lattice site, first it suppresses the 
possibly resident boson of the same type, since the site can carry multiple 
bosons only if they are characterized by different values {A, /u}- Then the 
new boson is assigned a momentum that equals the sample momentum 
of the particle. If WA is the stochastic variable describing the site boson 
momentum, a re-initialization of such variable occurs whenever a new boson 
pair is created according to the rule 

W£M <- “V = f («) 

One can interprete this rule as an exchange of momentum between the site 
and the particle (considering its sample momentum), through their respec¬ 
tive bosons. 

When the site is not visited or it senses a particle counter A = p (5 = 0), 
the momentum of the resident A/U-boson decays according to the rule 


(l) (t- 1 ) 

w \n = w \„ 


1 - 


(Sw 


( 0 ) \2 ' 
A/x > 




* = 1 , 2 ,..., 


(44) 


where £ is the lifetime of the resident boson, that is, the number of iterations 
spanned from when it was created (£ = 0 for a new boson). 

Rule (PHD applies separately to each boson type. Note that this rule may 
be interpreted as a discrete analogous of an exponential decay. Note also 
that, since A/r € Q and <5 ,£ £ N, also € Q. 


4.1.2 Quantum force on the particle 

Similarly to site bosons, when a new particle A/r-boson is created, it sup¬ 
presses the possibly carried boson of the same type, as the particle can carry 
multiple bosons only if they are characterized by different values {A, /u.}. 
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The particle boson momentum is assigned a value p'^ that equals that 
of the boson of the same type that was previously resident at the site, if 
there was one. If the momentum of the Xp boson possessed by the particle 
is represented by the new stochastic variable V\^ (the iteration index is 
omitted for simplicity), a re-initialization occurs whenever a new boson pair 
is created, according to the rule 

<- pZ = («) 

When the particle visits a site where new A/r-bosons are not created, the 
momentum of its A^-boson decays according to the rule 

P\,l = P&" ( X “ 2l) ’ k = l,2,... (46) 

where k is the lifetime of the particle boson, that is, the number of iterations 
spanned from when it was created (k = 0 for a new boson). Rule (1461) may 
be interpreted as a discrete analogous of a decay that is r\_/ l/Vk. 

Finally, the effective momentum propensity that regulates the particle 
motion is not the particle momentum p but its total momentum 

V T -=P-Y, V ^ ( 47 ) 

where the summation is taken over all bosons carried by the particle (the k 
index has been removed from the right-hand side terms for simplicity since, 
at a certain time step, each boson has its own lifetime). 

Simple as it is, the mechanism given by (f4hl) - (T471) is capable of accounting 
for self-interference, as it will be proven in next section for a number of 
increasingly complex situations. 

4.2 Ensembles of Particles 

In this section the predictions of QM are retrieved for ensembles of similarly 
prepared particles. The scenarios considered are (i) the two-slit or two- 
source scenario, (ii) the multiple-source scenario, (iii) the particle-in-a-ring 
scenario, and ’iv) the particle-in-a-box scenario. 

4.2.1 Two slits 

With the mechanism illustrated in Sect. sn derive now the probability 
function P(£, r) for the two-source scenario, that is, two lattice sources 
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separated by 5 sites, labelled 1 (source at 5/2) and 2 (source at —5/2), 
respectively. The probability that a particle is emitted from source 1 is Pi, 
while P 2 = 1 — Pi is the probability of emission from source 2. 

In this scenario, any site register CX can take only the values £ — 5/2 
(particle coming from source 1) and £ + 5/2 (particle coming from source 2). 
Let me label these two alternative events “1” and “2” in the following. Four 
types of Xfi events are possible, namely, “11”, “12”, “21”, and “22”. The 
probabilities of these events are easily calculated as Pn := P(, P 12 := P 1 P 2 , 
P 21 := P 1 P 2 , P 22 := P ‘2 ■ The events “11” and “22” do not create any boson 
(the particle counter is equal to the site register). Consider now the bosons 
“ 12 ”. 

The boson momentum at a site {£, r} is a stochastic variable whose 
expected value is calculated in Appendix I A. 71 The final result is 


E\Wi i2 ] 



(48) 


a result remarkable due to the emergence of a trigonometric function from 
the discrete-valued model proposed. 

On the other hand, the expected value of the particle boson momentum 
is calculated in Appendix IA.8l as 


E/P 12 ] — \fp\P 2 ■ 


sin(7r5^) 

ird 


(49) 


Again, the result is remarkable since a square-root function emerges from 
the discrete-valued model proposed. 

Finally, the expected value of the particle total momentum is obtained 
from (l47|) as 


e[v t ]= p -2^pXp 2 


sin(7r5d+) 

7 t5 


(50) 


since bosons “12” and “21” give the same contribution to the combined 
momentum (while the events “11” and “22” do not generate bosons). 

The average particle motion is thus described for large r’s by a couple of 
equations involving £ r := lim^oo E[X T ] and p T := linx^oo E[P T \, namely, 


£r+l £t + Pt 


(51) 


and 


p T = p-2 v / P, P 2 • 


sin 


7T<5 — 
t 


7 t5 


(52) 
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Using the properties E[X T ] « E[C T \ = tE[P t \, the steady-state prob¬ 
ability density f(p T ) can be calculated given the probability density of p, 
f(p) = 1/2. Using (1521) and the rule f(p T )dp T = f(p)dp, obtain 

f(pr) = = ^(1 + 2\/ PiP 2 cos(Tt5p T )). (53) 

2 dp T 2 

The probability density function P(£, r) still obeys the approximated rule 
(|35j) . that is 

P{(,T) = m=(/T) . (54) 

T 

From (1531) obtain finally 

1 + ‘1\J P\ P -2 cos (nd-j 

mr) = - 2r V T) ■ (55) 

This result is in perfect agreement with QM predictions, derived in Ap¬ 
pendix IA.6.31 from Schrodinger’s equation, under the equivalence rule P ~ 
IT| 2 derived in Sect. [3] and after transformation to lattice units. 


4.2.2 Multiple slits 

Similar considerations apply for scenarios with N s > 2 sources, where Pi is 
the probability of the i-th source (X^i Pi = 1) and $ij is the distance in 
lattice units between sources i and j, for i j. In these cases, there are 
N S (N S — 1) possible types of bosons, each of which ultimately contributes 
to a term y / PiPj cos | ^ to the momentum pdf. Since contribution of 
boson ij is the same as that of boson ji, they can be lumped to yield 


f(Pr) 


1 + X^i X^+i ^y/PjPj cos (^\djj\Pr) 


(56) 


1+Xi=i Xj=i+i ivm™ ( n \dij\ 'j 

P&r) = --- - - ^(57) 

Again, this result is in perfect agreement with QM predictions, derived in 
Appendix I A. 6. 41 from Schrodinger’s equation, after transformation to lattice 
units. 
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4.2.3 Particle in a Ring 


Also the well-known “lattice in a ring” scenario can be simulated with the 
presented model, since the particle has no interactions with the environment. 
Now the useful lattice does not extend infinitely in both directions but only 
to values £ E [0, •£ — 1], where £ is the ring circumference in lattice units. Due 
to ring periodicity, it is now necessary to distinguish the stochastic variable 
X T , the site reached after t time steps, from the particle spatial counter C T 
defined by (fl3l) . In fact (assume for simplicity Xq = 0), X T = C T mod £. 

The register at a site £ can take the obvious value £, but also the 
value £ + £, if the particle has reached the site after having crossed the entire 
ring once, the value £ + 2£, if the ring has been crossed twice, etc. Also 
values £ — £,£ — 2£, ... are possible, if the ring is crossed in the opposite 
direction. Summarizing, C% can take infinitely many distinct values £ ± n£. 
Consequently, all possible particle paths interefere, with a path difference 
that is always a multiple of £. In the proposed model this situation is 
equivalent to a multiple-source preparation, with infinitely many equally- 
probable and equally-spaced sources separated by a distance i in lattice 
units. 

In this scenario, the probability density function of the particle momen¬ 
tum is f(p') = 5{p' — p), that is, an ensemble of particles having the same 
momentum p is considered. The steady-state probability density of the total 
momentum p T is to be calculated. Equation (1521) is replaced by 


N s N s 


Pt=P - 


lim 

N s —yoo 


EE 

i=1 j—i-(-1 


2 sin (7r|j — i\£q) 
N s 7r| j - i\£ 


(58) 


The limit at the right-hand side of (j58p reads 


(1 

£ 

{ 1 ~^ + 

_^2 


(59) 


where |_-J denotes the “floor” function. Using this result and after manipu¬ 
lation equation (1551) is equivalent to 



p£ 

Y 


(60) 


where the operator [•] denotes the rounding function. In other words, for 

a well-defined particle momentum p, the steady-state momentum p T can 

2n 

only take one of the discrete values —, that for n = [^-]. This result 
coincides precisely with QM predictions derived in Appendix IA.6.51 after 
transformation to lattice units. 
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4.2.4 Particle in a Box 

Although in the “particle in a box” scenario the particle is not free but 
submitted to external forces as soon as it hits the box boundaries, this case 
is easily treated with purely geometric considerations. 

The spatial domain is limited to £ E [0,-£], where IX = L is the phys¬ 
ical width of the box. When the particle hits one of the boundaries, it is 
“reflected” due to the infinite potential that is assumed outside the box. In 
the proposed model, a reflexion (infinite external force) means two things: 
change of sign of (i) the momentum propensity, of the momentum of all 
bosons possessed, and of the spatial counter of the particle, and (ii) reset of 
its time counter to zero. Due to the latter change, the particle can interfere 
with its previous paths through creation of bosons. In fact, the register £$ 
at a site £ can take the obvious value £ (assume for simplicity £o = 0) if 
the site is visited before any reflexion, the value —i —{£ — £) if it is visited 
after one reflexion, the value 21 + £ if it is visited after two reflexions, etc. 
Consequently, the path difference is always a multiple of 21. 

In the proposed model this situation is equivalent to a multiple-source 
preparation, with infinitely many equally-probable and equally-spaced sources 
separated by a distance 2£ in lattice units. From considerations that are sim¬ 
ilar to those presented for the particle in a ring case, it is found that the 
steady-state probability density of momentum p T has peaks for the values j , 
with n = \p£\. This result coincides precisely with QM predictions derived 
in Appendix A.6.6, after transformation to lattice units. 

4.3 Numerical results 

This section is aimed at reproducing numerically the self-interference scenar¬ 
ios analysed in the previous sections. The model described in the previous 
sections is in principle implemented by the algorithm shown in Table [2j With 
respect to the algorithm of Table [1} more instructions are now needed to 
represent the creation of pairs of bosons, their initialization, and the decay 
of their momentum. 

Numerical results of Fig. [3] have been obtained with a modified version 
of the algorithm above, aimed at fastening calculations. In fact, the inter¬ 
ference pattern starts to build when the frequency of A/r events approaches 
the theoretical probabilities P \ M introduced in Sect. 14.21 That requires that 
a sufficiently large number of lattice sites are visited by a sufficiently large 
number of particles. This “lattice training” process requires in turn that 
a large number, say Npi, of particles are emitted. When the number of 
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Algorithm 2 Pseudocode used for the simulations of Fig. [3j 


9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 


for particle 1 to Np do 


> Initialization: 


£ random value bewteen —5 and +<5 
r <r~ 0 

pO 4— random value beween -1 and +1 
A «- 0 

for all possible bosons b do 
p(b) <r~ 0 

k\b) <- 0 

end for 

for iteration 1 to Nt do 


> Particle dynamics: 


T i — T + 1 

p p0 — J2p(b) 

v random value +1, 0 or -1 with prob. given by (jlOl) 

c«- £ + ^ 

A <— A + v 
q ^ X/t 

> Momentum decay for particle bosons, see (j46]): 
for all possible bosons b do 
k(b) <- k(b) + 1 
p(b) ■(— p{b)(l - 1/2 /k(b)) 

end for 

> Momentum decay for lattice bosons, see (1441) : 
for all sites s do 

for all possible bosons b do 
l(s, b ) <r- l(s , b) + 1 

w(s,b) <— w(s,b ) * (1 — (w0(s,b)/l(s,b)) 2 ) 

end for 
end for 

> Boson creation: 

site s 4— {£, r} 
if p(s) — A = —6 then 
boson b «- “12” 
else if p(s) — A = 6 then 
boson b <r- “21” 


else 

no boson b 

end if 


> Boson momentum reinitialization: 


p(b) <— w(s, b) 
k(b) <- 0 
w(s, b) <— q 
wQ(s, b ) 8q 
l(s, b) •(— 0 
end for 

(A,p(s)} ^ {p,(s), A} 
v(£) <- + 1 
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Figure 2: Values of w \2 in lattice of size Np after Np = 100 particle emis¬ 
sions, vs. respective steady-state values. 

particle emissions required for the statistical build up of the intereference 
pattern is considered, it becomes clear that complete simulations would re¬ 
quire extensive computing times. 

To fasten computations, I have considered the “lattice training” already 
completed. In practice, the values /r found by a particle are not acessed in 
the site registers but randomly attributed with probabilities that equal the 
steady-state ones. In the two-slit scenario, ^ equals q\ with probability Pi, 
and q 2 := (£ + <5/2)/r with probability P 2 . Similarly, particle bosons are 
directly initialized with their expected steady-state momenta, see (TIH1) . In 
fact, the convergence of the process (11071) requires much less iterations than 
the training of the lattice, and it is completed well before Np± emissions. 
The latter statement can be observed in Fig. [2] that shows pairs of values 
W 12 and gisinc( 2 ( 7 i) in a lattice after Np = 100 and for Np = 300 (no gap 
between two successive emissions). 

Figure [3] shows the frequency calculated for the two-slit scenario 
(5 = 2) after a time Np = 300 for different values of Np, with the lattice 
already trained as specified above. As the the number of particles emitted 
in the ensemble increases, a frequency distribution builds up. For large Np, 
the frequency clearly tends to the a priori probability P(^,Np) given by 

flSSD- 

Two more scenarios are presented in Fig. |4j The former is a case of two 
sources at distance 5 = 2 with non-equal probabilities Pi = 0.9, P 2 = 0.1. 
The second scenario is that of three equally-probable sources located at £ = 
{—1,0,1}. The frequency values have been obtained with the accelerated 
algorithm discussed above. In the first scenario, two types of bosons are 
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Figure 3: Frequency of arrival of particles emitted at £ = ±1 as a function 
of £ after Np = 300. From top-left to bottom-right, Np = 500, 5000, 10000, 
and 50000, respectively. 


created, both having creation probability P 1 P 2 and distance 6. In the second 
scenario, there are two types of bosons with <5 = 2 and four types of bosons 
with (5 = 4, all of them with probability equal to 1/9. In both cases, the 
figure shows a precise agreement between large r’s prediction of the model 
with the theoretical probability densities. 

5 Discussion and Future Work 

The proposed approach has been proven capable of describing trajectories 
of individual particles in an ensemble of similarly-prepared particles in a 
simple and realistic way. Simple and realistic means that the ontology of 
the proposed model includes real particles, a real discrete spacetime lattice, 
both capable of storing and exchanging a few pieces of information, and 
arithmetic operations. In fact, all quantities modeled are represented either 
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x 1 o 4 ^ x 1 o 4 


Figure 4: Number of particle arrivals for: two sources at £ = {±1} with 
respective probability 0.1 and 0.9 (left), three equally-probable sources at 
£ = {—1,0,1}; both scenarios with Np = 10000, Np = 5000; red curves: 
theoretical probability densities. 


by integers or rational numbers, and operations on them are products and 
sums. 

The predictions of the model have been shown to tend to the predictions 
of QM in the continuum limit for free particles and, more remarkably, also in 
the case of quantum interference. In my opinion, the latter evidence makes 
the model a successful candidate to provide both qualitative and quantitative 
explanation for several quantum phenomena. However, other QM aspects 
have still to be added to the model. 

The results of this paper concern free particles only. However, the pro¬ 
posed model seems naturally capable to integrate also external forces into 
the picture. Each interaction of the particle with its sourrounding is indeed 
expected to modify its intrinsic properties. External forces shall be treated 
as permanent variations of the particle intrinsic momentum propensity p. 
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due to new types of bosons carrying a specific momentum that amounts to 
that of the particle. The particle counters shall be reset similarly to the 
particle-in-a-box scenario. Relativistic Newton’s second law shall be con¬ 
sidered, to prevent that p becomes larger than unity under the action of 
persistent forces. 

Extension to two- and three-dimensional spaces seems also natural. A 
set of three momentum propensities, p x , p y , and p z shall be introduced, 

f ulfilling the condition that the total energy e = - ! Vx+ ^" ' v ~ < 1. This 
condition implies that p^. + p 2 y + p\ < 1, thus fixing a constraint to the 
probability densities of the three propensities. 
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A Appendices 

A.l Frequency and matter waves 


In some interpretations of quantum phenomena, a particle is associated 
with a matter wave, whose frequency is proportional to its energy via 
the Planck constant. In the proposed model, the frequency is retrieved 
as the reciprocal of the average return time to any position of the lattice. 
To see that, define the probability mass function P(n ) as the probability 
that a particle returns at an arbitrary position for the first time after 
a time 2 n. For example, P{ 1) = 4ac = b 2 / 2, P( 2) = 2a 2 c 2 + 2 ab 2 c = 
5/86 , etc. The general expression for P{n) is 


P(n) = 2b 2n 



n 


+ E 


1 + 4(n 

4 & 



2 b 2n 


fn 4 13 

V3 ~~ 9 + ~9~ 



( 61 ) 
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Now, define the average return time as 


Using the results 

OO 

J2 kr 

k=1 

one can find that 


T r := E[n] = 


i. r 

- (T^T)2- 


E”=i nP ( n ) 

i P{n) 


(62) 


k=\ 


2 k = r{l + r) 
(1 — r) 3 ’ 


(63) 


n =1 


" 3 (1 — 6 2 ) 2 9 V(l-6 2 ) 2 


E p < 

O 

E rf ( 


2 b 2 


-in? 


9 V(! — & 2 /4) 2 


- 1 


2 6 2 (1 + 6 2 ) 8 b 2 26 6 2 /4 

U) = r—--TTTTTT + 


n=l 


3 (1 — 6 2 ) 3 9 (1 — 6 2 ) 2 9 (1 — 6 2 /4) 2 ’ 


(64) 

(65) 


and consequently r r as a function of b. Let me now introduce the energy 
with the substitution b = 1 — e. After some tedious but straightforward 
manipulations of (!Ml) - (l65l) . find the frequency / = l/r r as 


f = e(2 - e)(—1 - e)(3 - e)(e 4 - 4e 3 + 5e 2 - 2e + 1) 

1 (e 2 — 2e — l)(5e 4 — 20e 3 + 29e 2 — 18e + 6) 1 J 

that is the relationship sought. It is easy to verify (Fig. [5]) that for 
e = 0, / = 0, while for e = 1, also / = 1. Moreover, for small values of 
e, the relationship (l66h is approximated by 


/ = e, 


(67) 


which is precisely the de Broglie relation in lattice units. 


A.2 Alternative derivations of (I18j) 

Consider (ITUl) as a binomial distribution f(k;n,q ) with k = £ + r, 
n = 2 t, q = (1 + p)/ 2. For n large enough an approximation of p 
is a normal distribution with mean p = nq = t + rp and variance 
er 2 = nq{\ — q) = br. Recalling that k — p = £ — pr, obtain (fT8l) . 
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Figure 5: Matter wave frequency as calculated from (1661) and as f = e. 


Yet a third possible method would start from expressing the recursive 
equation (fT5j) as 


pAO 


p r _l(C + l)-Pr-l(e-l) , 
Pr-l(£) = ~P -j- + 

, Pr-1 (Z + 1) - 2 pr-l (£) +PT-AZ-1) 
+ e ' -^-• 


( 68 ) 


The latter difference equation has a continuum limit described by the 
differential equation 

dp dp d 2 p 

d t = ~ P dti + 6 W 



which is a convective-diffusion equation with e playing the role of the 
diffusivity and p of the convection velocity. Note, however, that in (1181) 
the correct result for the diffusivity is b = 1 — e and not e as it would 
be predicted by ([69]) . 


A.3 Special Relativity 

Equations (jl6l) or (118(1 are invariant with respect to Lorentz transforma¬ 
tions. To see that, let me take a particle having momentum propensity 
p and a reference frame that is moving at velocity V with respect to 
the fundamental lattice. In this refercne frame, the fundamental lattice 
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dimensions (primed quantities) are deformed as follows: 


, , = X{l-qP)y/rrp = T(l-qP)y/T=W 

{l-p/3) 2 ’ ‘ (1 -p/3 ) 2 


( 70 ) 


where <7 = 7 and f5 := V/c, so that the relationship X'/T' = X/T = c 
still holds. 

Consequently, 




(1 -P/3) 2 
(l-/3 2 )(l-g/3) 


(£- 0t), 


/ 


T 


(1 -pP? 

(1 — /3 2 )(1 — qP) 


m (7i) 


so that 


x' = i'X' 


y(x — Vt), t' = t'T' = 7(4 



(72) 


with 7 = 1 / 1 — /3 2 , in agreement with Lorentz transformations. 

In both reference frames, the probability of having the particle at a 
site of the fundamental lattice (£, r) must be the same. In the primed 
reference frame, (fl8l) is calculated as 


ptfS) 


1 ( (£'- p'r'f\ 

V 2 ^ 6 V eXP l 26V )■ 


(73) 


Now, let me assume that the momentum propensity in the moving 
reference frame is 


/ P~P 

P = 


1 - Pp 


(74) 


in agreement with relativistic velocity addition formula. Consequently, 
since b' = (1 — p' 2 )/ 2 , one easily verifies that 


6 ' 


, 1 -P 2 

(l-p/3) 2 ' 


(75) 


With these relationships, it is easy verified that 


and 


- P'r' = (£ - pr) 


1 -pP 

1 ~qP 


b'r' = 6t 


(76) 

(77) 
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Using (|76l) and (1771) . it is easily verified that the quantity p( t') given 
by (1731) approximates p(£, r) given by (1181) when q « p, that is, for large 
r’s. In reality, all the formulae above are valid for large r’s, since, for 
instance, I do not impose that £' or t' are integers. 

In our laboratory system (primed reference frame moving with an un¬ 
known velocity V with respect to the fundamental lattice), we prepare 
an experiment with a momentum propensity that we label as p' and we 
observe a probabilty mass function at a point in the spacetime that we 
label as x', t'. Although our assumptions on these values is incorrect, 
we measure the correct probability mass function using the p formula. 

If (3 = p, that is, in the reference frame of the particle, we obtain the 
following results: X' = jX, V = 7 T, £' = £ — pr, r' = r — p£, p' = 0, 
b' = 1/2. Thus the particle itself “sees” a broader lattice ( X' > X, 
T' >T). For /3 = p —>• ±1, the lattice becomes infinitely large. 

A.4 Derivation of (12711 

The number of paths leading to a certain site is given by 



(78) 


where n a , and n c are the number of moves with v = 1,0,—1, re¬ 
spectively. Clearly, n a + rib + n c = t and, in order to reach exactly the 
site in question, n a — n c = £. With these constraints, (fTSl) becomes 



(79) 


The range of n c derives from the fact that all the binomial arguments 
are positive integers. Therefore, n c +l; > 0, r— n c —£ > 0, r—2n c —^ > 0. 
Consequently, max(0, —^) < n c < . 

The probability of each path is given by a na b nb c nc . Since b 2 = 4ac, 
this probabilty is also calculated as 2 nf >a ria+Tli 7 2 c nc+Tli 7 2 . From (fT6l) . the 
probability of all N p paths is 



2 r 

T + C. 


T+e, t -4 

a 2 c 2 



a n a +nf,/2 C n c +ni,/2 

(80) 


30 






For each value of n c there will be a number N p (n c ) of paths with the 
same cumulated energy. The value taken by <S| equals the sum n a +n c = 
2 n c + £ since only these moves contribute to the cumulated energy. The 
relative probability of these N p (n c ) paths over the totality of N p paths 
is thus (1271) . 

For the sake of completeness, one might define a particle stochastic 
variable as well, 

n n 

Sn := X] l V «'l = E V «'’ ( 81 ) 

n'=riQ n'=no 

whose pmf follows the binomial distribution with r trials and e = a + c 
probability of success, i.e. 

(j> T (a) := Pr(5 T = a) = ^^\e a b T ~ a . (82) 

with support a = {0,..., r} (e is the energy propensity, not the Neper 
number). 


A.5 Derivation of (12311 

From the definition ([22 ]) . 


pi 



2 2t+1 


2 2t+1 B(t + £ + 1,t — £ + 1) 


where B(-) is here the Beta function, 


B(v,w) := f 1 z v - x ■ (1 - z^dz. 

Jo 

From the properties of such function, it follows that 


( 2r \ (r + £)!(r — £)! 

U + £ j (2r + l)! 


(83) 


(84) 


(85) 


that is, (1231) . 
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A. 6 Probability densities for the scenarios considered 
A.6.1 Single source 

The wavefunction for a free particle is 

*(x,t) = —L [ e l{kx -'^v(kS))dk, (86) 

v2vr J 

where the wavenumber k is related to the momentum of the particle 
and 

<p(k, 0) = -j=f *(®', 0 )e- ikx 'dx'. (87) 

For a single perfectly localized source at x = 0, (1B71) reads 


<p(k,o) = 

and consequently (|86j) is rewritten as 

X 


X 

y/2ir 


\£ f(x,t ) = 


that, integrated, yields 




f 

Hk 2 

/ exp 

i(kx - 2m^ 


dk 


^(x, t) = 


X 





exp 


2i 


mx 
1 2 ht 


y/2n\ iht 

The probability density is easily calculated as 

|T , ..o mX 2 mX 2 


2irht 


ht 


( 88 ) 


(89) 


(90) 


(91) 


A.6.2 Equally-probable sources 

The solution of the Schrodinger equation for the case where the two 
sources are equally probable is based on the linear superimposition of 
the two waveforms, 




’J'l (x,t) + T 2 (x,t) 

V2 


(92) 
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where both \&i(a;, t) and ^ 2 (x,t) are obtained from (1901) by replacing 
the term x 2 (that was valid for a source at x = 0) with a term (x — ^X) 2 
and (x + | X ) 2 , respectively. Thus, 


/ m 1 

mix — ix) 2 

+ exp 

1 

l O|C< 

+ 

1_ 

\ 

exp 

1 2 ht 

2 ht 

i 





(9 


T(x,t) = 


The probability density is given by 

2 mX 2 


l *(M)r = 


ht 


cos 


5i - 5, 


(94) 


where S\ and S 2 are the two independent action values, that is, the 

‘ItttyiSXx 

phases of the two exponentials in (1931) . As S\ — S 2 = -—-, obtain 


l*0M)| 2 = 


mX 2 

ht 


1 + cos 27r 


m5Xx\ 

~ht ) 


(95) 


and thus an interference term arises with respect to (1251) . due to the 
presence of two possible sources. The interference is related to the phase 
difference between the two waveforms. 


A.6.3 Non equally-probable sources 

In the case the probability of the two sources is different, say, Pi 7^ 
P2, equation (l92jl is replaced by T(x’,f) = (,x,t). 

Consequently, (1951) is replaced by 

l>Kx.<)| 2 = ^( 1 + 2/fWco S (2,=A£)) (96) 

If Pi = P2 = 1/2, (|95l) is retrieved. 


A.6.4 Multiple sources 

In the case of N s sources, (f95l) is generalized as follows. 
mX 2 


l*(M)r = 


ht 


N s N a 

i+E E 2 v qs x COS | 27T 
i=l j=i -\-1 


m\Xi — Xj\x 
ht 
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where x % is the locations of the z-th source and '}2n=\ -Pj = 1 • If the 
sources are equally probable and equally spaced by a distance 5X , (|95l) 
is particularized as 


l*M )| 2 


mX 2 

ht 


1 + 


Ns -1 

E : 

j= 1 


N s -j 

N s 


cos 



mj5Xx\ 
ht ) 


(98) 


A.6.5 Particle in a Ring 


The “particle in a ring” scenario can be described by assuming a single 
plane wave of a well-defined wavenumber and momentum. Equation 
(186(1 is still valid with ip(k, 0) = 5(k), where 5 is here the Dirac delta 
function (assume positive direction of momentum). As the term e lut 
does not affect the amplitude of the wavefunction, it can be neglected 
to write 

(®) = ~^e ikx . (99) 

>/2n 

However, periodicity imposes that 'h(x) = 'h(x ± 2nR), where R is the 
ring radius, or e lkx = e * fc (z+ 27r ^), which implies that e lk2nR = 1. This 
condition is satisfied for 2irkR = 27rn, where n E Z, or for 


Tl 

kn = ( 100 ) 

in agreement with the somehow standard derivation that considers an¬ 
gular momentum and polar coordinates. In terms of momentum p = Hk 
(here the dimensional momentum, not the momentum propensity as in 
the paper body), (HOOD reads 

nh 

Pn = T , 



where L = 2ixR is the ring circumference. 


A.6.6 Particle in a Box 

The wavefunction has the same form ()102l) of the “particle in a ring”, 
however, with opposite amplitudes for the forward and backward waves, 

T(x) = -^=(,°, ikx - e~ ikx ), (102) 
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Var. 

Support 

Exp. value 

Variance 

pmf 

V n 

v € {0,±1} 

p 

i -p“ 

2 

{a, b , c} 

Xn 

£ € {£o — T, . . . , £o + t} 

Co + pr 

l~p\ 

2 1 

M£) 

C n 

A e {-t, ... ,t} 

pr 

l ~P\ 

2 

Pn(£l 0 + A) 


£ 

£ 

0 

Deter. 

of 

{0,1} 

M£) 

MO(i-MO) 

(MO, i- 

S n 

cr e {0,... ,r} 

1 +P T 

2 

i?(l-i±2)r 

M 0 ") 

sf 

a e {|£l, l£l + 2, • • •, |£| + 2 2 } 

ew-r 

2t—1 

(J29D 

<t>T 


satisfying the boundary condition \k(0) = 0. The second boundary con¬ 
dition 'I'(L) = 0 imposes that e lkL = e~ lkL or e 2lkL = 1. The latter 
condition is satisfied for 2kL = 2irn, where n € Z, or for 

, n7r 

kn = (103) 


In terms of dimensional momentum, p n 


nh 

2L 


holds. 


A.6.7 Step Potential 
A.7 Derivation of (14811 


I will consider 12-bosons and omit all boson (“12”) and site (£,r) in¬ 
dices for the sake of simplicity. The boson momentum at a site {£, r} 
is a stochastic variable depending on its lifetime, another stochastic 
variable, according to the update rule 


W V- 


w(°\ 



Pr = r' 

Pr = 1 - r' 


l <- 


0 , 

£ + 1 


Pr = r' 

Pr = 1 — r' 


(105) 


where the probability r' that a new site boson is created is the product 
of three terms: the frequency with which particles in the ensemble 
are emitted, p £ (the probability that the site is visited), and P \2 (the 
probability of the event 12). 


Mm 
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The quantity w ® is defined by rule (1431) . Clearly, it is constant after 
that the site is visited for the first time and 


w® = qi 


(106) 

holds, with qi := ^^. Thus the quantities are also constant 
and can be calculated by repeated application of rule 


w (t} = qi JJ 

3 =1 


1 - 


A 2 


(107) 


The probability that the site carries a boson with lifetime £ and mo¬ 
mentum equals the product r'(l — r') 1 Consequently, the expected 
value of W is 


E[W] 




(108) 


Since the values of r' are generally very small, (11081) is rewritten for 
large r’s as 

E[W] = q\ JJ ^1 - = 9iSinc(5gi). (109) 


However, for large r’s q\ ~ which is true also for <72 
Therefore, for both bosons “12” and “21”, 


E[W] 


sin(7r<5 
7 t5 




holds. 




( 110 ) 


A.8 Derivation of (14911 

The probability that the particle carries a boson “12” equals the proba¬ 
bility that it carries a new boson, or a boson that is one iteration-old, or 
a boson that is two iterations-old, etc., all these events being mutually 
exclusive. Such probability is therefore calculated as 

T—1 

Pr[boson“12” = True] = Pi 2 +Pi 2 ^-Pi 2 )+Pi 2 {l-Pi 2 ?+- ■ ■ = Pu^il-Pn) k 

k =0 

(111) 
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However, we consider large r’s probabilities to retrieve P(£, r), so let r 
in (11111) tend to infinity. 

Boson lifetime affects the momentum exchange, as per (1461) . Thus the 
mean value of the boson momentum is calculated as 


E[P 12 ] = Pl 2 ]T(l-Pi 2 ) 


VS? 


( 112 ) 


k=0 


From the definition (1461) . 


Ik) k 

— = TT 

( 0 ) 11 

Pl2 1=1 


p V2 tt 21—1 1 • 3 • 5 • 7 • 


21 


2-4-6' 


(2fc)! 

(kl) 2 4 k ' 


(113) 


After some manipulations, the latter formula turns out to be formally 
equivalent to 


(fc) 

Pl2 _ (_ U k 
(0) ^ ’ 
P\2 


- 1/2 

k 


(114) 


as it can be verified by inspection. Therefore, (11121) is rewritten as 

E{Va] - u 2 f;(i-Pi2) i (-i)‘P/ 2 ) = p.2(i-(i-pi 2 ))- i/2 = vT? 

(115) 


(o) 

PV2 


k=0 


after having recognized in the right-hand side the binomial series 


(i-x)« = f; 


k=0 


a 


(-i)' 


(116) 


multiplied by P\ 2 - 

The quantity ’jiy] is a stochastic variable on its own, since it takes a 
different value (that of the resident boson) for each site visited when 
a new pair of bosons is created. From (l45l) and (l48|) this stochastic 
variable tends to be distributed as 


(0) _ sin(TrJ^) 

Pl 2 ~ ttS ' 

With this result, equation (1 115 D is rewritten as 

sin(7rh—) 


E[P 12] = • 


tt5 


(117) 


(118) 
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